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A  COMPARATIVE  STUDY  OF  LINEAR  APRAY  SYNTHESIS 
TECHNIQUE  USING  A  PERSONAL  COMPUTER 


INTRODUCTION 

Procedures  for  synthesizing  the  radiation  patterns  of  linear  arrays 
based  on  the  specification  of  its  si  delobe  structure  are  well  establish¬ 
ed.  One  of  these  is  the  technique  originally  proposed  by  Do1plr  and 
it  provides  for  a  uniform  si  delobe  level.  Another  technique,  although 
developed  for  line  sources,  is  due  to  Taylor^  and  can  be  adopted  for 
linear  arrays.  The  latter  is  popular  due  to  its  synthesized  aperture 
distributions  which  are  more  readily  realizable.  A  discrete  version  of 
the  Taylor  synthesis  procedure  is  discussed  by  Shelton-^. 

Both  Dolph-Chebyshev  and  Taylor  synthesis  techniques,  fundamental¬ 
ly,  rely  on  manipulation  of  the  zeros  of  the  linear  array  pattern 
function.  The  aperture  distribution  for  the  desired  pattern  function 
usually  requires  lengthy  computations.  In  case  of  Dolph-Chebyshev 
synthesis,  this  problem  has  been  addressed  by  several  authors4-^ 
over  the  past  few  decades.  The  Taylor  synthesis,  in  effect,  uses  a 
discrete  Fourier  Transform  technique  (called  WoodwardlO  synthesis) 
to  obtain  the  aperture  distribution.  These  procedures  do  not  have 
much  in  common;  as  a  matter  of  fact,  in  case  of  an  endfire  Chebyshev 
array,  expression  for  the  element  excitations  are  quite  different  from 
that  for  a  broadside  Chebyshev  array.  However,  the  knowledge  of  the 
pattern  null  locations  in  the  above  synthesis  procedures  can  be  used 
to  develop  a  simple  expression  that  is  suitable  for  all  cases.  The 
expression  is  readily  developed  based  on  the  convolution  synthesis 
procedure  discussed  by  Laxpati 11  for  planar  arrays. 

With  several  alternate  expressions  being  available  for  the  aperture 
distribution  of  Dolph-Chebyshev  and  Taylor  syntheses,  it  is  desirable 
to  undertake  a  study  to  make  some  recommendations  as  to  the  suitability 
of  these  expressions  in  numerical  computation.  Due  to  the  increasing 
use  of  personal  computers  by  antenna  engineers,  it  is  felt  that  an  in¬ 
vestigation  of  this  nature  should  be  confined  to  the  computation  using 
such  small  computers.  Thus,  in  this  paper,  we  present  the  results  of  a 
comparative  study  of  various  linear  array  synthesis  techniques.  In  the 
next  section,  after  a  brief  discussion  of  the  three  basic  techniques 
for  evaluation  of  Chebyshev  coefficients,  we  discuss  the  accuracy  and 
computation  times  associated  with  these  techniques.  The  following 
section  presents  the  results  of  the  study  involving  two  different 
techniques  (one  due  to  Shelton^  and  the  other  using  the  convolution 
procedure)  for  Taylor  synthesis.  In  the  last  section,  some  general 
observations  about  the  investigation  and  on  the  results  are  offered. 

DOLPH-CHEBYSHEV  SYNTHESIS 

Following  Dolph's  paper  on  Chebyshev  synthesis,  Barbiere^, 
Van  Der  Maas5,  Salzev^,  and  Brown?, 8  reported  on  alternative  means  of 
evaluating  aperture  distribution  for  Chebyshev  arrays.  Although  they 
are  not  the  same,  the  expressions  by  Barbiere,  Salzev  and  Brown  are 
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similar  in  that  they  express  the  current  in  an  element  in  terms  of  a 
finite  series  of  terms  involving  ratios  of  factorial  functions  and  with 
alternating  sign.  The  expression  by  El  1 i ott 12  is  representative  of 
this  group  and  is  the  one  used  in  this  work  and  is  reproduced  below. 
We  shall  call  this  the  classical  expression.  Also,  although  our  results 
are  valid  for  odd  or  even  number  of  elements,  for  simplicity,  we  will 
present  examples  of  odd  number  of  elements.  Thus,  all  linear  arrays 
discussed  in  the  following  have  (2N+1)  elements;  the  element  numbering 
scheme  is  shown  in  figure  1,  where  the  elements  are  assumed  to  have  a 
symmetric  excitation  leading  to  the  broadside  radiation. 

Classical  Technique: 


y  rnN-p-^ 

L  ("1/  N+p 
p=n 


rfN+p+i)  («0) 

r(N-p+l)r(p+n+l)r(p-n+l) 


n=0,l,2,... ,N. 

where  T2n(u0)=R  and  SLL  =  20  log  R.  Here,  SLL  is  the  desired  sidelobe 
level  in  dB,  T2n(x)  is  the  Chebyshev  polynomial  of  degree  2N  and  r(x) 
is  the  Gamma  function. 

In  contrast,  the  expression  given  by  Van  Der  Maas  involves  terms 
of  the  same  sign  inside  the  summation.  Bresler9  reformulated  the 
expression  into  a  recursive  form  using  nested  products.  This,  we  feel 
is  a  distinctly  different  form  of  representation  of  the  coefficients. 
Thus,  we  use  this  representation  (called  Nested  Product  Technique)  in 
our  comparison.  This  expression  (in  our  notation)  is  shown  below. 

Nested  Product  Technique: 

iN-n  =  2N  a  NP(n,fm,a).  n=0,l,2,  (2) 


where  NP(n,fm,a)  =  l  an-m  n  f j ;  fn  =  1. 
m=l  j=m 

.  .  m(2N-2n+m) 

and  fm  =  (n-m)(n+i-m)  '» 


also  a  =  1 


The  third  technique  is  based  on  the  convolution  of  three  element 
canonical  arrays*!.  These  canonical  arrays  have  outer  element  excita¬ 
tions  of  unity,  whereas  the  center  element  excitation  cj;  for  j =1 ,2 , 
...N  is  chosen  such  that  the  jth  canonical  array  has  a  pattern  null 
at  the  location  of  the  jth  symmetric  zero  pair  of  the  Chebyshev 


polynomial.  These  arrays  are  then  convolved  to  generate  the  large 
array.  Convolution  Technique: 


Woj  -  cos  2N  ^  »  j=l  ,2,. . . ,N# 

Where  w0j  are  the  zeros  of  the  Chebyshev  polynomial  T2n(w). 

Cj  =  -2cos  u0j ; 

u0j  =  2arccos  (w0j/u0). 

And  the  aperture  distribution 
N 

I(x)  =  I  In6(x-nd)  =  *fN  (3) 

n=-N 

where  fj  =  6(x-d)+CjS(x)  +  6(x+d). 

Using  these  three  expressions  (equations  (1),  (2)  and  (3)),  comput¬ 
er  programs  NESTED,  CHEB  and  CONCHEB,  respectively,  were  written  to 
implement  the  Chebyshev  synthesis.  Different  versions  of  the  program 
suitable  for  implementation  on  different  machines  were  written.  These 
were  two  personal  computers  used  in  the  numerical  phase;  one  is  an  8- 
bit  Radio  Shack  TRS-80  Model  II  which  has  available  an  interpretive 
BASIC  language.  The  other  computer  is  a  16-bit  NEC  Advanced  Personal 
Computer  with  BASIC  and  FORTRAN  IV  compilers.  Also,  in  order  to 
ascertain  the  numerical  accuracy,  some  of  the  programs  were  run  on  a  32- 
bit  mainframe  computer  (Texas  Instrument's  Advanced  Scientific  Computer 
at  the  Naval  Research  Laboratory)  using  double  precision  (REAL*8) 
arithmetic. 

The  computation  was  carried  out  for  several  different  array  sizes 
ranging  from  15  to  99  elements;  although,  in  principle,  there  is  no 
limit  to  the  size  of  arrays  that  may  be  synthesized.  Furthermore,  all 
designs  specified  a  si  delobe  level  of  30  dB. 

Figure  2  shows  the  run  time,  under  FORTRAN,  for  the  three  afore¬ 
mentioned  Chebyshev  synthesis  programs  versus  number  of  elements.  The 
CHEB  program  was  the  slowest;  but  more  importantly,  the  program  failed 
to  converge  to  the  correct  element  excitations  beyond  30  elements. 
Over  21  elements  the  accuracy  of  the  excitation  was  only  to  2  digits. 
When  the  program  was  run  using  double  precision  arithmetic  it  still 
failed  to  converge  above  31  elements.  This  indicates  that  the  classical 
technique  inherently  has  a  limitation  as  to  the  largest  size  of  array 
that  may  be  synthesized. 

The  convolution  synthesis  program,  CONCHEB,  although  much  faster 
than  CHEB,  certainly  cannot  complete  with  NESTED  program  in  speed. 
Also,  beyond  61  elements,  the  CONCHEB  program  failed  to  converge.  The 
current  version  of  the  program  convolves  three-element  arrays  using 
zeros  of  the  Chebyshev  polynomial  in  an  alternating  sequence;  i.e.,  the 
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Figure  2  -  Run  Times  for  Chebyshev  Programs  - 
FORTRAN  IV  on  NEC-APC 
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sequence  in  j  is  1,N,2,N-1,3, — .  No  attempt  was  made  to  modify  this 
convolution  process  to  improve  the  accuracy;  past  experience  with  the 
convolution  process  indicates  that  some  improvement  may  be  possible. 
However,  with  reference  to  figure  2,  it  is  obvious  that  the  NESTED 
program  is  the  most  efficient  one. 

The  results  of  the  element  excitations  indicate  that  it  is  an 
extremely  stable  algorithm;  provides  a  very  good  accuracy  in  single 
precision  (six  digit  accuracy);  and  of  course,  it  is  very  fast.  This 
program,  NESTED,  was  translated  into  BASIC  and  run  on  the  TRS-8D,  Model 
II  computer.  The  execution  time  ranged  from  two  seconds  for  21  element 
array  to  36  seconds  for  a  99  element  array.  Although,  the  execution 
times  in  BASIC  are  about  15  to  20  times  longer  than  that  in  FORTRAN, 
they  are  not  significantly  long  to  be  of  any  major  consequence.  The 
NESTED  program  was  also  run  using  double  precision  (16  significant 
digits)  on  the  mainframe  computer.  The  total  execution  time  for  all 
10  different  arrays  was  less  than  0.3  seconds! 

Our  experience  with  synthesis  of  various  Chebyshev  arrays  using 
these  three  different  techniques  clearly  demonstrates  that  the  most 
important  consideration  on  small  computers  is  not  the  speed  of  execution 
but  the  accuracy  of  the  final  result.  In  this  sense  as  well,  the  nested 
product  algorithm  proposed  by  Rresler9  is  the  winner. 

TAYLOR  SYNTHESIS 

Synthesis  procedure  proposed  by  Taylor2  applies  to  a  continuous 
aperture.  In  practice,  this  procedure  is  used  for  discrete  aperture 
(arrays)  by  properly  discretizing  the  continuous  distribution.  Shelton3 
presented  a  synthesis  procedure  for  discrete  aperture  distribution  for 
Taylor  type  sidelobe  structure.  He  expressed  the  pattern  function  in 
the  form  of  a  product  function  of  zeros  and  then  carried  out  the  synthe¬ 
sis  exactly  analogous  to  that  by  Taylor;  that  is,  to  use  the  Woodward 
synthesis  technique.  In  particular,  for  a  2N+1  element  array,  all  2N 
zeros  are  explicitly  specified  in  the  pattern  function.  Thus,  analogous 
to  the  Chebyshev  synthesis,  this  synthesis  is  amenable  to  the  convolu¬ 
tion  procedure.  In  view  of  this,  in  the  case  of  Taylor  synthesis,  we 
compare  the  two  techniques;  one  proposed  by  Shelton  and  the  other 
being  the  convolution  synthesis.  Before  presenting  and  discussing  the 
results  of  the  investigation,  the  pertinent  expressions  for  the  two 
syntheses  are  given  below.  Once  again,  we  will  limit  out  discussion 
to  arrays  with  odd  (2N+1)  number  of  elements. 

Discrete  Taylor  (Shelton3)  Technique: 

2tt7T  y/  A2+(n-l/2)2 

uon  =  (2N+1)  i  ♦  n=l>2,  —  ,n-l 

“V  A2+(n-l/2)2 


2irn 

=  (2N+1) 


,  n=n,  — ,N. 


(4) 


I 

where  A  =  it  cosh-1 (R);  n  is  equal  to  the  number  of  near-in  zeros  that 
are  moved  in  order  to  achieve  the  desired  sidelobe  ratio  R  (or  equiva¬ 
lently  the  number  of  near-in  sidelobes  that  are  required  at  the  speci¬ 
fied  level).  The  element  excitations  are 


n-l  2mpir 

Ip  =  1+2  l  am  cos  2N+1  .  p=0,— ,N  .  (5) 

m=l 

where 

.  2nm  . 

am  =  E  l  2N+1J  '* 


N 

E(u)  =  n 


n=l 


(cosu-cosuon) 

(l-cosuon) 


For  the  case  of  the  convolution  synthesis  procedure,  once  the 
symmetric  zero  pairs  are  established,  the  excitation  of  the  center 
element  of  a  three  element  canonical  array  is  readily  determined.  The 
procedure  and  expressions  are  analogous  to  the  case  of  Chebyshev 
convolution  synthesis.  They  are 

zeros  are  +  u0j;  j=l,2,-~,N 

where  u0j  are  defined  through  equation  (4),  and  the  excitation 

Cj  =  -2cosu0j. 

The  synthesis  of  the  large  array  is  carried  out  using  the  convolu¬ 
tion  of  three  element  arrays,  chosen  in  the  same  alternating  zero 
sequence  as  indicated  for  the  Chebyshev  array. 

Based  on  these  two  procedures,  computer  codes  STAYL  and  CONTAYL, 
respectively,  were  developed  in  FORTRAN  using  single  precision  arithme¬ 
tic  Run  time  associated  with  these  codes  for  rT  =  6  and  the  sidelobe  level 
of  30  dB  for  various  number  of  elements  from  15  to  99  were  recorded  and 
are  shown  in  figure  3. 

The  program  CONTAYL  failed  to  converge,  once  again,  for  arrays 
with  more  than  71  elements  and  provided  only  two  to  three  digit  accuracy 
between  31  and  61  elements.  These  results  are  similar  to  the  Chebyshev 
convolution  synthesis.  Even  the  run  time  data  is  very  close. 

The  computation  time  associated  with  STAYL  has  an  interesting 
behavior  with  increasing  number  of  elements;  it  is  almost  linear.  This 
is  to  be  expected,  since  the  number  of  computations  to  be  carried  out 
for  each  element  is  determined  by  ”n  and  not  (2N+1).  The  corresponding 
growth  for  CONTAYL  is  exponential.  Thus,  for  small  number  of  elements 
CONTAYL  may  save  some  computation  time  but  will  suffer  in  accuracy  as 
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TIME  IN  SECONDS 


11  21  31  41  51  61  71  81  91  101 


NO.  OF  ELEMENTS 

Figure  3  -  Run  Times  for  Taylor  Programs  - 
FORTRAN  IV  on  NEC-APC 
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the  number  of  elements  increases.  A  check  of  STAYL  program  using  double 
precision  arithmetic  on  the  mainframe  computer  indicates  that  it  has 
five  to  six  digit  accuracy  in  single  precision  on  a  small  computer. 

It  should  be  noted  that  the  STAYL  program  code  was  developed  by 
Shelton  for  the  HP-41C,  a  pocket  calculator.  On  this  calculator,  one 
has  10  significant  digit  capacity  and  thus  the  results  obtained  are 
more  accurate  than  with  a  single  precision  FORTRAN  Code.  But,  as  one 
would  expect,  the  HP-41C  is  very  slow;  it  took  approximately  5  minutes 
to  synthesize  a  31  element  array. 

STAYL  Code  was  also  run  on  NEC-APC  using  CBASIC,  a  compiler  BASIC. 
In  CBASIC,  the  computation  times  were  significantly  higher,  ranging 
from  30  seconds  for  a  15  element  array  to  217  seconds  for  a  99  element 
array.  However,  the  computation  was  carried  out  to  14  significant 
figures. 

Once  again,  as  with  Chebyshev  synthesis,  we  find  the  overriding 
consideration  in  Taylor  synthesis  is  not  the  computation  time,  but  the 
accuracy  of  the  results.  In  this  sense,  Shelton's  procedure  is  most 
efficient. 

CONCLUSIONS 


As  is  often  the  case  with  engineering  investigations,  the  most 
significant  results  presented  in  this  paper  are  not  what  we  were  looking 
for  when  we  began  the  project.  We  were  originally  interested  in 
evaluating  computer  run  times  for  the  various  programs.  However,  two 
points  soon  became  apparent  —  first,  most  of  the  programs  run  fast 
enough,  even  on  small  machines,  so  that  run  time  is  not  a  major  concern, 
and  second,  only  two  of  the  programs  give  adequate  precision  for  the 
range  of  array  size  that  was  investigated.  It  is  concluded  that 
Bresler's  nested  product  algorithm  gives  excellent  results  in  terms  of 
speed  and  precision,  and  also  that  Shelton's  discretized  procedure 
allows  precise  Taylor  synthesis  for  all  sizes  of  arrays.  Finally,  it 
is  noted  that  the  programs  are  very  brief;  the  FORTRAN  computer  codes 
for  all  five  programs  are  included  in  the  appendix  and  the  codes  in 
BASIC  are  also  available  from  the  authors. 
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APPENDIX 


In  this  appendix  the  FORTRAN  IV  computer  codes  for  NEC-APC  with 
supersoft  FORTRAN  compiler  are  listed.  As  noted  in  the  main  body  of 
the  report,  the  programs  are  brief;  there  are  a  number  of  "comment" 
statements  in  the  listing  and  thus  are  easy  to  follow.  No  sample  inputs 
or  outputs  are  included. 

Programs  listed  in  the  following  are  NESTED,  CHEB,  CONCHEB,  STAYL 
and  CONTAYL. 


71 


001 

C 

PROGRAM  NESTED 

002 

C  REVISED  01/28/84 

002 

C  CHEBYSHEV  ARRAY  ALGORITHM  USING  NESTED  PRODUCTS  FORMULATION 

003 

C  INPUT  IS  M  =  NO.  OF  ELEMENTS;  SLL  =  SIDELOBE  LEVEL  IN  DB. 

004 

REAL  1(100),  NP,  CC 100) 

005 

WRITE  (1,100) 

006 

100 

FORMAT  (’  ENTER  DATA:  M , SLL ’ ) 

007 

READ  (1,200)  M,SLL 

008 

200 

FORMAT  (  I0,F0.0) 

009 

N  =  M  /  2 

010 

TEST  =  (-1 )*#M 

011 

IF(TEST.GT.O)  GO  TO  10 

012 

N  *  (M  -  1)  /  2 

013 

10 

R  =  10.  **  (  SLL  /  20.  ) 

014 

ARCOSH  =  ALOG  (  R  +  SQRT  (  R»»2.  -  1  )  ) 

015 

A  =  ARCOSH  /  (M-1 ) 

016 

ALPHA  =  (  TANH  (A)  )  »«  2. 

017 

I(N+1 )  =  1.0 

018 

I(N)  =  (M-1)  *  ALPHA 

019 

DO  30  K  =  2,  N 

020 

NP  =  1.0 

021 

DO  20  J  =  1,  K-1 

022 

FN  =  J  *  (M-1-2»K+J) 

023 

FD  =  (K-J)  •  (K+1-J) 

024 

F  =  FN/FD 

025 

NP  =  NP  «  ALPHA  •  F  +  1. 

026 

20 

CONTINUE 

027 

I(N+1-K)  =  (  M-1  )  •  ALPHA  »  NP 

028 

30 

CONTINUE 

029 

DO  40  L  *  1,  N+1 

030 

C(N+L)  =  I(L) 

031 

C(N+2-L)  =  I(L) 

032 

40 

CONTINUE 

033 

WRITE  (4,  50) 

034 

50 

FORMAT  (  '  CURRENTS’) 

035 

DO  60  L  =  1 ,  M 

036 

WRITE  (4,70)  L,  C(L) 

037 

60 

CONTINUE 

038 

70 

FORMAT  (10X,  12,  10X,  F10.6) 

039 

STOP 

040 

END 

12 


001 

002 
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006 

007 

008 

009 
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011 

012 

013 

014 

015 

016 

017 

018 

019 

020 

021 

022 

023 

024 

025 

026 

027 

028 

029 

030 

031 

032 

033 

034 

035 

036 

037 

038 

039 

040 

041 

043 

044 

045 

046 

047 

049 

050 

051 

052 


C  PROGRAM  CHEB 

C  REVISED  01/28/84 

C  BASED  ON  A  CLASSIC  METHOD  OF  COMPUTATION  OF  CHEBYSHEV 
C  EXCITATION  VOLTAGES. 

C  REFERENCE  ANTENNA  THEORY  AND  DESIGN;  ELLIOTT. 

C  ODD  NUMBER  OF  ELEMENTS  ONLY 

REAL  CC(100),C(100) ,CRNT( 100) 

WRITE  (1,100) 

100  FORMAT  ('  ENTER  DATA:  N,  SLL' ) 

READ  (1,200)  N.SLL 
200  FORMAT  (IO,FO.O) 

M  =  N-1 

MM  =  M/2 

NN  =  (N+1)/2 

PI  =  3-1415927 

R  =  10.  »•  (SLL/20.0) 

U  =  COSH  ( RCOSH  (R)/  FLOAT  (M)  ) 

DO  10  I  =  1 ,NN 
II  =  I  -1 
C(I)  =  0.0 
DO  20  J  =  I,NN 
JJ  =  J  -  1 
A  =  FLOAT  (NN  +  JJ) 

GA  =  GAMALN  (NN  +  JJ) 

GB  s  GAMALN  (NN  -  JJ) 

GE  =  GAMALN  (J  -  II) 

GD  =  GAMALN  (J  +  II) 

UP  s  U  **  (2»JJ) 

SIGN  =  (-1)  «•  (NN-J) 

TL  =  EXP  (GA  -  GB  -GE  -  GD) 

TN  =  UP*SIGN*(2. *NN-1 . )/(2.*A) 

T  =  TL  •  TN 
C(I)  =  C(I)  +  T 
20  CONTINUE 

10  CONTINUE 

DO  12  J  =  1 ,NN 
CC(J)  =  C(J)/  C(NN) 

12  CONTINUE 

DO  13  J  =  1 ,NN 
CRNT  (NN-1+J)  =  CC(J) 

CRNT  (NN+1-J)  s  CC(J) 


13 

CONTINUE 

WRITE  (4,30) 

30 

FORMAT  (» 

DO  14  I  =  1,N 

CURRENTS') 

WRITE  (4,300) 

I,  CRNT(I) 

14 

CONTINUE 

300 

FORMAT  (10X,I2 

,15X,F12.8) 

STOP 

END 


13 


051 

C* 

052 

C 

INVERSE  HYPERBOLIC  FUNCTION 

053 

C* 

054 

FUNCTION  RCOSH  (R) 

055 

RCOSH  =  ALOG(R  +  SQRT(R*R  -  1.0)) 

056 

RETURN 

057 

END 

058 

C» 

059 

C 

HYPERBOLIC  FUNCTION 

060 

C» 

061 

FUNCTION  COSH  (R) 

062 

Y  =  EXP  (R) 

063 

COSH  =  (Y  +  (1.0/Y)  )  /2. 

064 

RETURN 

065 

END 

066 

C» 

067 

C 

GAMALN  FUNCTION 

068 

c» 

069 

FUNCTION  GAMALN  (K) 

070 

GAMALN  r  0.0 

071 

IF  (K  .EQ.  0)  RETURN 

072 

FACT  =  0.0 

073 

TPL  =  0.91893853 

074 

AL  =  K 

075 

10 

IF  (AL  .GE.  10.0)  GO  TO  20 

076 

FACT  =  FACT  +  ALOG  (AL) 

077 

AL  s  AL  +  1.0 

078 

GO  TO  10 

079 

20 

TERM  =  (AL  -  0.5)  •  ALOG(AL)  -  AL  +  TPL 

080 

1  +  1. 0/(12. «AL)  -  1.0/(360.0  »  AL**3)  • 

081 

2  AL**5)  -  1.0/  (1680.  «  AL**7) 

082 

GAMALN  =  TERM  -  FACT 

083 

RETURN 

084 

END 

14 


1. 0/(1260.  • 
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PROGRAM  CONCHEB 


001 


002 

c»»< 

003 

C* 

LINEAR  ARRAY  SYNTHESIS  USING  CONVOLUTION  METHOD. 

004 

C» 

CHEBYSHEV  SIDELOBE  DESIGN 

005 

C» 

ODD  NUMBER  OF  ELEMENTS. 

006 

>•«•« 

007 

REAL  PSI(IOO) 

008 

REAL  COOOMAaoOMKIOOMEOOOMSOOO.COHVOOO) 

009 

DATA  A1/100*1 ./,  A3/100*1 ./ ,  C0NV/100»0./ ,  AA/100*0./ 

010 

IHH 

Oil 

c* 

N  =  NUMBER  OF  ELEMENTS  IN  THE  ARRAY.  MUST  BE  ODD! I 

012 

c» 

SLL  =  SIDE  LOBE  LEVEL  IN  DBS. 

013 

c»«< 

»«»#» 

014 

PI  =  3.1415297 

015 

WRITE  (1,100) 

016 

100 

FORMAT  (•  ENTER  DATA:  N,SLL«) 

017 

READ  (1,200)  N,SLL 

018 

200 

FORMAT  (IOjFO.O) 

019 

M  =  N-1 

020 

NN  =  (N+1)/2 

021 

MM  =  (N-1)/2 

022 

MD  =  (MM/2)  +  1 

023 

20 

CALL  CHEBX ( PI , M , NN , SLL , MM , PSI ) 

024 

025 

C* 

CREATE  THREE  ELEMENT  ARRAYS 

026 

IHHI 

027 

30 

DO  40  -I  =  1 ,MM 

028 

A2(I)  =  -2.  *  COS(PSKI)) 

029 

40 

CONTINUE 

030 

CM* 

IHIH 

031 

C» 

REPEATED  CONVOLUTION  OF  3-ELEMENT  ARRAYS 

032 

C»«H 

033 

CONV(1)  =  A1(1) 

034 

C0NV(2)  =  A2( 1 ) 

035 

C0NV(3)  =  A3(1) 

036 

L  =  1 

037 

K  =  5 

038 

LX  =  0 

039 

50 

LL  =  NN-L 

040 

LX  =  LX+1 

041 

60 

L  =  LL 

042 

C(1)  =  A1(L) 

043 

C(2)  =  A2(L) 

044 

C ( 3)  =  A3(L) 

045 

DO  70  I  =  1,3 

046 

DO  70  J  =  I,K 

047 

JJ  =  J-I+1 

048 

AA( J)  s  AA( J)  +  CONV( JJ) *C(I) 

049 

70 

CONTINUE 

050 

DO  80  I  =  1,K 

051 

CONV(I)  =  AA(I) 

052 

AA(I)  =  0.0 

053 

80 

CONTINUE 

054 

K  =  K+2 

055 

IF  (L.EQ.MD)  GO  TO  90 

056 

LL  *  NN+1-L 

057 

LSUM  x  L+LX 

059 

IF  (LSUM.EQ.NN)  GO  TO  60 

060 

GO  TO  50 

061 

90 

CONTINUE 

062 

WRITE  (4,600) 

063 

600 

FORMAT  (  •  CURRENTS'/) 

064 

DO  120  I  x  1,N 

065 

WRITE  (4,700)  I,CONV(I) 

066 

120 

CONTINUE 

067 

700 

FORMAT  (1  OX, 12,1  OX, FI 0.6) 

068 

STOP 

069 

END 

070 

FUNCTION  COSH(R) 

071 

Y  =  EXP(R) 

072 

COSH  =  (Y  +  ( 1 .0/Y) )/2. 

073 

RETURN 

074 

END 

075 

c*«< 

IHH* 

076 

C* 

INVERSE  HYPERBOLIC  COSINE  FUNCTION 

077 

»*•»«* 

078 

FUNCTION  RCOSH(R) 

079 

RCOSH  =  ALOG(R  +  SQRT(R*R  -  1.0)) 

080 

RETURN 

081 

END 

082 

c*»* 

»*••** 

083 

C» 

CHEBYSHEV  ZEROS 

084 

imi* 

085 

SUBROUTINE  CHEBX ( P I , M , NN , SLL , MM , PSI ) 

086 

REAL  X(50),PSI(100) 

087 

R  =  10.0  ««(SLL/20.) 

088 

B  =  COSH(RCOSH(R)/M) 

089 

DO  10  I  a  1 ,NN 

090 

J  =  1-1 

091 

X(I)  a  C0S(PI*(2.*J  +  1 . )/(2.*M) ) 

092 

10 

CONTINUE 

093 

DO  20  J  a  1  ,MM 

094 

II  a  NN-1+J 

095 

JJ  a  NN-J 

096 

Y  a  X(J)  /  B 

097 

PSI(II)  a  2.»ATAN(SQRT(1-Y»Y)/Y) 

098 

PSI(JJ)  a  PSI(II) 

099 

20 

CONTINUE 

100 

RETURN 

101 

END 

16 


001 

C 

PROGRAM  STAYL 

002 

C  REVISED  01/28/84 

003 

C  THIS  PROGRAM  COMPUTES  ELEMENT  EXCITATIONS  FOR  TAYLOR 

004 

C  TYPE  SIDELOBES  USING  SYNTHESIS  EXPRESSIONS  OF  SHELTON 

005 

DIMENSION  2(100)  ,AM( 100) ,EN( 100) ,EX( 100) 

006 

WRITE  (1,100) 

007 

100 

FORMAT  ('  ENTER  N , NBAR , SIDELOBE  LEVEL  FOR  TAYLOR’ 

008 

1  SYNTHESIS'  ) 

009 

READ  (1,200)  N,  NBAR,  SLL 

010 

200 

FORMAT  (2I0,F0.0) 

Oil 

WRITE  (4,300)  N,  NBAR,  SLL 

012 

300 

FORMAT  (’  TAYLOR  SYNTHESIS  -  SHELTON’  /  ’  N=',I5, 

013 

1  2X, ’NBAR=',I5,2X,'SIDEL0BE=’,  F5.2) 

014 

AL2  =  0.30102999566398 

015 

ALE  =  0.43429448190325 

016 

PI  =  3.14159265358979 

017 

M  =  (N-1)/2  +  0.1 

018 

IE  =  1 

019 

IF  (N  .EQ.  (2*M+1 ) )  IE  *  0 

020 

A  =  (SLL  +  20.0*AL2)  /  (20.0*PI*ALE) 

021 

XN  =  FLOAT  (N) 

022 

XN12  =  FLOAT  (NBAR)  -  0.5 

023 

N1  =  NBAR  -  1 

024 

ALPHA  =  SQRT  (A«A  +  XN12  »  XN12) 

025 

DO  1  1=1, N1 

026 

X1 12  =  FLOAT  (I)  -  0.5 

027 

BETA  =  SQRT  (A*A  +  XI12  *  XI12) 

028 

Z(I)  =  ( (2.0*PI/XN)/ALPHA)  *  FLOAT  (NBAR)  •  BETA 

029 

1 

CONTINUE 

030 

DO  2  I=NBAR ,M, 1 

031 

Z(I)  =  (  FLOAT  (I)  »  2.0*PI)/XN 

032 

2 

CONTINUE 

033 

EO  =  1.0 

034 

DO  3  1=1,  M 

035 

3 

EO  =  EO  *  (1.0-COS(Z(I))) 

036 

DO  5  1=1, N1 

037 

AM(I)  =  1.0 

038 

DELTA  =  (2.0*PI  *  FLOAT(I) )/XN 

039 

IF  (IE  .EQ.  1)  AM(I)  =  COS  (DELTA/2.0) 

040 

DO  4  J=1,M 

041 

4 

AM(I)  =  AM(I)  *  (COS (DELTA)  -  COS(Z(J))) 

042 

AM(I)  =  AM(I)/EO 

043 

5 

CONTINUE 

044 

DO  6  1=1, M+1 

045 

XI  =2*1-2 

046 

IF  (IE  .EQ.  1)  XI  =  XI  +  1 

047 

EN(I)  =  0.0 

048 

DO  7  J=1 ,N1 

049 

XJ  =  FLOAT  (J) 

050 

EN(I)  =  AM( J)  «  COS( (PI*XI*XJ)/XN)  +  EN(I) 

051 

7 

CONTINUE 

052 

EN(I)  =  2.0  •  EN(I)  +  1.0 

053 

6 

CONTINUE 

054 

DO  50  K=1 ,M+1 

055 

L  =  N+1-K 

056 

EX(K)  =  EN ( M+2-K ) /EN ( M+ 1 ) 

057 

50 

EX(L)  =  EX(K) 

058 

WRITE  (4,301) 

059 

WRITE  (4,55)  (I,EX(I),  1=1, N) 

060 

301 

FORMAT  ('  ELEM.  NO.',  3X,  ‘EXCITATION’) 

061 

55 

FORMAT  (5X,I2,5X,F14.7) 

062 

STOP 

063 

END 

001 

C 

PROGRAM  CONTAYL 

002 

c»»< 

>•••• 

003 

C» 

LINEAR  ARRAY  SYNTHESIS  USING  CONVOLUTION  METHOD. 

004 

C* 

TAYLOR  SIDELOBE  DESIGN 

005 

C* 

ODD  NUMBER  OF  ELEMENTS. 

006 

nil* 

007 

REAL  PSI(IOO) 

008 

REAL  C(100),AA(100),A1(100),A2(100),A3O00),C0NV(100) 

009 

DATA  A 1 / 1 00*1 ./,  A3/100*1 ./,  CONV/100*0./,  AA/100*0./ 

010 

C««) 

IH»I 

011 

c» 

N  =  NUMBER  OF  ELEMENTS  IN  THE  ARRAY.  MUST  BE  ODD! ! 

012 

c» 

SLL  a  SIDE  LOBE  LEVEL  IN  DBS. 

013 

IHM 

014 

PI  =  3.1415297 

015 

WRITE  (1,100) 

016 

100 

FORMAT  ('  ENTER  DATA:  N, NBAR,SLL’ ) 

017 

READ  (1,200)  N,NBAR,SLL 

018 

200 

FORMAT  (2I0.F0.0) 

019 

M  =  N-1 

020 

NN  =  (N+1)/2 

021 

MM  a  (N-1)/2 

022 

MD  =  (MM/2)  +  1 

023 

10 

CALL  TA YLX( MM , SLL , PI , PSI , N , NN , M , NBAR ) 

024 

imii 

025 

c* 

CREATE  THREE  ELEMENT  ARRAYS 

026 

027 

30 

DO  40  I  a  1 ,MM 

028 

A2(I)  a  -2.  «  COS(PSKI)) 

029 

40 

CONTINUE 

030 

c»»< 

mm 

031 

c« 

REPEATED  CONVOLUTION  OF  3-ELEMENT  ARRAYS 

032 

>•**•» 

033 

CONV(1)  a  A1(1) 

034 

C0NV(2)  a  A2(  1 ) 

035 

CONV( 3)  =  A3(1) 

036 

L  a  1 

037 

K  =  5 

038 

LX  a  0 

039 

50 

LL  a  NN-L 

040 

LX  a  LX+1 

041 

60 

L  a  LL 

042 

C(1)  =  A1(L) 

043 

C(2)  a  A2(L) 

044 

C(3)  =  A3(L) 

045 

DO  70  I  a  1,3 

046 

DO  70  J  =  I,K 

047 

JJ  a  J-I+1 

048 

AA( J)  a  AA( J)  +  CONV(JJ)*C(I) 

049 

70 

CONTINUE 

050 

DO  80  I  a  1,K 

051 

CONV(I)  8  AA(I) 

052 

AA(I)  a  0.0 

053 

80 

CONTINUE 

054 

K  a  K+2 

055 

IF  (L.EQ.MD)  GO  TO  90 

056 

LL  =  NN+1-L 

057 

LSUM  =  L+LX 

058 

IF  (LSUM.EQ.NN)  GO  TO  60 

059 

GO  TO  50 

060 

90 

CONTINUE 

061 

C 

WRITE  (4,401) 

062 

C  401  FORMAT  (  '  PSI  ZEROS'/) 

063 

C 

DO  110  I  =  1,M 

064 

C 

WRITE  (4,500)  I,PSI(I) 

065 

C  110  CONTINUE 

066 

500 

FORMAT  (10X,I2,10X,F10.6) 

067 

WRITE  (4,600) 

068 

600 

FORMAT  (  '  CURRENTS'/) 

069 

DO  120  I  =  1 ,N 

070 

WRITE  (4,700)  I,CONV(I) 

071 

120 

CONTINUE 

072 

700 

FORMAT  (10X,I2,10X,F10.6) 

073 

STOP 

074 

END 

075 

FUNCTION  COSH(R) 

076 

Y  =  EXP(R) 

077 

COSH  =  (Y  +  (1.0/Y))/2. 

078 

RETURN 

079 

END 

080 

081 

c» 

INVERSE  HYPERBOLIC  COSINE  FUNCTION 

082 

083 

FUNCTION  RCOSH(R) 

084 

RCOSH  =  ALOG(R  +  SQRT(R*R  -  1.0)) 

085 

RETURN 

086 

END 

087 

c»«< 

088 

c» 

COMPUTATION  OF  TAYLOR  ZEROS 

089 

c««< 

090 

SUBROUTINE  TAYLX ( MM , SLL , P I , ZEROS , N , NN , M , NBAR ) 

091 

REAL  ZERO(50),  ZEROS(IOO),  MEMA , MEMB 

092 

NBAR1  =  NBAR-1 

093 

A  =  (SLL  +  6.0202)/27.2875 

094 

095 

C»  COMPUTE  ZEROS  FROM  1  TO  NBAR 

096 

c««< 

097 

DO  10  I  =  1 ,NBAR1 

098 

RI  =  I 

099 

MEMA  =  (A*A)  +  ( (RI-.5) **2 . ) 

100 

MEMB  =  (A*A)  +  ( (NBAR-.5)**2. ) 

101 

ZERO(I)=( ( (2.*PI)*NBAR)/N)*( (SQRT(MEMA) )/(SQRT(MEMB) 

102 

10 

CONTINUE 

103 

c»« 

104 

C» 

COMPUTE  ZEROS  FROM  NBAR  TO  M 

105 

c« 

106 

DO  20  I  =  NBAR, MM 

107 

RI  =  I 

108 

ZERO(I)  =  (2*PI*RI)/N 

20 


109 

20 

CONTINUE 

110 

C 

WRITE  (4,100) 

111 

C 

100 

FORMAT  (  *  ZEROS 

112 

c 

DO  30  I  =  1 ,MM 

113 

c 

WRITE  (4,200)  I,ZERO(I) 

114 

c 

30 

CONTINUE 

115 

c 

200 

FORMAT  ( 10X, 12 , 10X,F10.6) 

116 

DO  40  J  =  1 ,NN 

117 

ZEROS(NN-UJ)  =  ZERO(J) 

118 

ZEROS(NN-J)  =  ZERO(J) 

119 

40 

CONTINUE 

120 

RETURN 

121 

END 
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